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The behavior of electronic states of one dimensional correlated disordered systems 

OO ; 

CN \ which are modelled by a tight binding Hamiltonian is studied analytically using the 
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. invariant measure method. The approach of Bovier is generalized to include the 

possibility of different site energies and nearest neighbor hopping integrals inside the 
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correlated sites or the cluster. The process is further elaborated by applying to the 
symmetric random trimer model which contains in it many hitherto known models of 
this category. An alternative mathematical definition of the exceptional energy (Eg) 
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Ch , from the invariant measure density along with physical arguments substantiating 
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it is presented. Furthermore, the procedure for obtaining exceptional energies is 
outlined and applied to the symmetric random trimer model to derive conditions for 
obtaining doubly degenerate exceptional energies. The Lyapunov exponent ( , ~f(E)) or 
the inverse localization length of states around the exceptional energy is found to vary 
as ~ (E — Es) 2n in the leading order, n denotes the degeneracy of the exceptional 
energy. The density of states at the exceptional energies are calculated. We further 
propose that one dimensional correlated disordered systems can be mapped to a Lloyd 
model in which the width of the distribution of site energies is determined by the 
reflection coefficient of the cluster embedded in the lattice of the another constituent. 
The importance of our results is discussed. 
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1 Introduction 



One of the well established results in condensed matter physics is that all electronic 
eigenstates of a disordered one dimensional system are exponentially localized irre- 
spective of the strength of the disorder. The early work of Anderson |l] on uncor- 
related site diagonal disorder in the tight binding model (TBM) and of Mott and 
Twose form the basis of this result. Of course, the result of Anderson and Mott 
and Twose cannot be rigorously valid in one dimensional systems in which the dis- 
order is correlated. For example, in the context of a TBM it has been shown that 
correlated off-diagonal disorder || cannot localize the state at the band center. An- 
other example in this category is the model proposed by Dunlap, Kundu and Phillips 
(DKP) 0. The most well known example in the context of a TBH, however, is 
the random dimer model (RDM) ||. This is basically the offspring of the original 
DKP model. The generalization of the RDM requires the extension of the correlation 
beyond the nearest neighbor site and the introduction of different nearest neighbor 
hopping elements among the correlated sites or simply the cluster. Two good but 
simpler examples in this generalized category are the repulsive binary alloy (RBA) 
U and the symmetric random trimer model (SRTM) 0. The SRTM which is the 
further generalization of the RBA, yields two exceptional energies. The exceptional 
energy is the energy at which the reflection coefficient of the cluster embedded in the 
lattice of the other constituent vanishes. We further note that the number of excep- 
tional energy for a trimer embedded in a one dimensional lattice of another element 
cannot exceed two. Another interesting as well as important feature of the SRTM 
is that positions of these exceptional energies can be tuned by changing either the 
hopping element or the site energies of the cluster and these two energies merge as 
a limiting case. This has actually been shown by Giri, Datta and Kundu (GDK) 
[0] by conventional analysis of the reflection coefficient and by appropriate numerical 
simulations. Salient features of the SRTM can be verified by fabricating appropriate 
layered hetero junction and also by coupling quantum dots ||. Another potential area 
of application of this model is the field of organic conducting polymers ||. 

The most commonly used method for studying the electronic properties of one 
dimensional correlated systems is to analyze the reflection coefficient of the cluster 
in the neighborhood of exceptional energies. Since the system behaves like a weakly 
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disordered system in the neighborhood of these energies, a very good estimation, 
albeit not rigorous, of the Lyapunov exponent (inverse localization length) in these 
neighborhoods can be obtained from the reflection coefficient. On the other hand an 
estimation of the density of states (DOS) at an exceptional energy can in principle be 
obtained from the phase of the transmission coefficient of a single cluster through the 
Thoules' formula |]TD1 . To the best of our knowledge no such effort has been made in 
this direction. Hence, the DOS and the mean square displacement of the particle are 
calculated numerically Jll| to establish the presence of nonscattered states around 
these energies. 

Another way of looking at this state of affairs is the calculation of the Lyapunov 
exponent and the integrated density of states (IDOS) of the system from its invariant 
measure. In an attempt toward understanding the behavior of the IDOS of the one 



dimensional Anderson model in the weak disorder limit, Bovier and Klein [12, [L3| de- 
veloped a scheme for a perturbation expansion of the invariant measure of the model. 
From the modified perturbative expansion of the invariant measure, Bovier and Klein 
showed that at all energies E = 2 cosa7r with a rational, the IDOS of the Anderson 
model in the weak disordered limit has singularities. This was the extension of the 



results obtained previously by Kappus and Wegner |14] and Derrida and Gardner 
[|HJ. Bovier and Klein [|I^] further showed that for irrational a, their technique gives 
a unique invariant measure with finite coefficients to all orders of perturbation. This 
modified expansion has also proved to be a true asymptotic expansion of the invariant 
measure 



16fl . This scheme was later applied by Bovier |L7| to develop the pertur- 
bation series expansion of the invariant measure around the exceptional energies of 
the RDM. This enabled him to show that the Lyapunov exponent vanishes as e 2 in 
energy (e) in the neighborhood of the exceptional energies. Furthermore, the IDOS 
is found to vary as e within this energy width. This is basically the first rigorous 
calculation on the RDM confirming the results of Ref . || . 

The fundamental characteristic of the cluster correlated disordered systems is the 
presence of exceptional energies where disorder systems purportedly behave like per- 
fect systems. However, for these systems to play an important role in transport 
properties of the materials, there must be a finite DOS at these energies. Hence, to 
fully characterize these systems, we need rigorous analytical calculations of the Lya- 
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punov exponent around these energies and DOS at these energies. To the best of our 
knowledge, for these systems, the invariant measure technique is the only technique 
that can yield the DOS analytically without invoking any approximation. This is 
primarily the motive to apply this technique to the SRTM which encompasses many 
hitherto known examples. For this purpose, we generalize the approach of Bovier [T7J 
to include the possibility of different site energies and nearest neighbor hopping ele- 
ments in the cluster. We also give here an alternative mathematical definition of the 
exceptional energy from the invariant measure and derive from it in the case of the 
SRTM, an algebraic equation in energy involving relevant parameters of the system. 
We further show that the equation gives the correct prediction of the possibility of 
tuning of exceptional energies yielding in the limit a degenerate exceptional energy 
as noted by GDK 0. Finally we propose a mapping of these systems to an effective 
Lloyd model. Such a mapping will be useful in the study of transport properties of 
these systems. 

The organization of the paper is as follows. In the following section we generalize 
the approach of Bovier. We then develop equations for the invariant measure density, 
the Lyapunov exponent and the IDOS. In section 4 we discuss the method for obtain- 
ing exceptional energies . Section 5 is devoted to the perturbative calculation of the 
invariant measure density for the SRTM. In section 6 and section 7 we calculate the 
leading order behavior of the Lyapunov exponent and the DOS at exceptional ener- 
gies respectively. In section 8 we deal with the mapping of the aspect. We conclude 
the paper by highlighting the major contributions of the paper. 



2 Formalism 

We study the Hamiltonian 

H = a l a n + V n+l,n{ a li+l a n + ^Ifln+l) (1) 
n n 

on l 2 (Z) where a n (aD destroys (creates) a particle at the nth site. V n ^ n+ i is the 
tunneling matrix connecting the nth site to the (n + l)th site. {V^+i^} are taken to 
be real and positive, although this constraint is not necessary for the formalism. 
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The eigenvalue equation |18[ associated with H is 



e nC n + V n+ i tTl C n+ i + V^ )Tl _iC n _i — i?C n (2) 

We introduce z n = Vn ^~ lCn g R with R denoting the compactified real line RU {ooj. 
The recursion relation for z n is then 

— E — € n ■ = £E,enV n ,n-i( Z n) (3) 

Zn 

We further note in this connection that the eigenvalue equation for a one dimensional 
array of masses {m,} coupled to nearest neighbors by identical harmonic springs is 

(2 - rriiQ 2 )ui = Ui-i + u i+ i (4) 

where Ui is the displacement of the ith mass, in the vibration with frequency 
Q. All spring constants are taken to be unity without any loss of generality. Now 
introducing a variable, z n = -, — c G R, we obtain from (ffl) 

z n+l — + C n : (5) 

where e n = ^- and V n ,n-\ = (^n^n-i) -1 ^ 2 - Hence, the behavior of one dimensional 
array of masses coupled by harmonic springs is mathematically equivalent to the 
one dimensional quantum motion of a particle in a TBM [TJL |20]]. The Lyapunov 
exponent, "f(E) and the IDOS, N(E) are related to the large n behavior of z n . If we 
define 

1 N z 

7(E) = lim — Y In — — (6) 

then 

7 (£) = Re 7 (£) (7) 

and 

iV(£) = -Im 7(E) (8) 

7T 

To understand the origin of equation (§) we consider a chain of N sites with fixed 
boundary conditions at both ends. In other words we set C_i = Cm = 0. To keep 
the argument simple, we assume that n G Z are real. From equation (3) we 

get zi — E — ei, a real quantity. So, none of the {z n } can be truly complex. However, 
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I T/ Z " }, n G Z can be negative real numbers which can be thought of as complex 

*n . n — 1 

numbers with the minimum phase, ir. So, the right hand side of equation (||) picks 
up a contribution whenever -jt^ — = ■^ a - is negative. 

^ V n ,n-1 Cn-1 ° 

Consider now the case when C n = for n < N. From equation (3) we obtain 

T/ 2 T/ 2 

n— l,n— 2 7-i n— 2,n— 3 / A \ 

z„_i = — ! — E — e„_ 2 ! (9) 

& ~ Crt-l Zn-2 

Values of E for which this equation (|9|) is satisfied are the eigenvalues of the system. 
When n = iV, we obtain eigenvalues of the system under study. Let us assume that Ei 
is the /-th eigenvalue of the system in ascending order. So, the eigen vector belonging 
to Ei has (I — 1) nodes. This in turn implies that v Zn - will also pick up (/ — 1) 
negative numbers giving Im-L J2n=i y^ El ^ = On the other hand, the number of 
sates upto E = E t , i.e IDOS (E{) is I. So, in the limit N — > oo, equation @ yields the 
IDOS (Ei). To understand further the behavior of equation @ for Ei < E < Ei + i, 
we note that Ei(m) > Ei(N) if m < N. For E in this limit equation (Q) will be 
satisfied for some m such that I < m < N and the /-th eigen value of the reduced 
system will be obtained. Hence, the number of modes and consequently the number 
of negative values of ( T/ z " ) will be preserved. We further note that when E = Ei + e 

Vn,n — 1 

and e — > 0, m will be close to N. When E = Ei + \ — e, m will be close to I and it will 
swing back to N for E = Ei + \. For further discussion on this see ref . pi | . 



The disorder in the Anderson model can arise from the disorder in diagonal ele- 
ments (e n ) of H, from the disorder in off-diagonal elements, V n n +i, n G Z of H or 
from both. In all these cases equation(3) defines a Markov chain in which states are 
characterized by the random variable z n , n G Z. When V njn+ i = V , n G Z and {e n } 
are i.i.d random variables, this Markov chain consists of persistent non-null states 



2!|. In other words, the chain is ergodic. Furstenberg's theorem ]23| then asserts 
that e n 7^ e, e = a constant, n G Z, there is a unique invariant measure (Ive{z) on R. 
This measures satisfies 

/ du E (z)f(z) = E / du E {z)f{E -e--) (10) 

for all bounded measurable functions, /. Here, E denotes the expectation with respect 
to the probability distribution of e. Furthermore, this measure is actually continuous 
and hence has support on R. So, the measure dv E [z) has a density, i.e., dv E (z) = 
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<f>E(z)dz, where 4>e(z) defines the density at z G R. Similarly, if e n = e, n G Z 
and {V^ jn _|_i} are i.i.d random variables, the resulting Markov chain is also ergodic 
except at E — e. Hence, in this case also a unique invariant measure exists on R ||13|| . 
From these established results we conclude that when both {e„} and {V„,n+i} ar e 
i.i.d random variables, the resulting Markov chain also consist of persistent non-null 
states. So, again a unique invariant measure will exist on R. 

To apply these results to correlated disordered systems some modifications are 
needed. In order to introduce the required modifications, we first describe the model. 
The model considered here is a random binary mixture of two types of clusters. Each 
cluster contains q > 2 elements. In the host cluster all elements are assumed to be 
same while in the guest cluster at least one element, if not all, should be distinct from 
the host element. All site energies and nearest neighbor hopping integrals in the host 
cluster are set to zero and unity respectively without any loss of generality. In the 
guest cluster site energies and nearest neighbor hopping integrals are allowed to be 
different. We further assume that the hopping integrals between the end sites of any 
two clusters are unity. So, diagonal and off-diagonal elements in our model are not 
totally random. Instead each category is required to satisfy q constraint relations : 

Cgm+i-l = [Q-l Pqm + Cq-l (1 _ Pqm)] e qm (H) 

and 

V qm +l,qm+l-l = 1 — [1 — (Vl t l-l p qm + X^— i (1 — P qm ))\(l — 5 q j)e qm (12) 

where 1 < I < q and qm G Z. The randomness in the model is, therefore, introduced 
through two Rademacher variables, e gm ,p qm G {0, 1}. These variables by construction 
are indeed i.i.d random variables. Rademacher variables, {p qm } are introduced to take 
into account the asymmetry of the guest cluster. Since a cluster can take only two 
possible orientations in the lattice, two possible values of {p qm }, namely zero and 
unity, occur with probability |. 

Since diagonal elements of H here determine the strength of hopping to and from 
the sites, the off-diagonal elements are not truly random. In essence this model 
is similar to the Anderson model of uncorrelated site disorder. However, due to 
constraints on diagonal and off-diagonal elements of H, {Z n } as such do not form a 
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Markov chain. To form the required Markov chain we need to consider the clusters as 
unit cells. In other words, we need to define a new random variable {X n } such that 

2*n+l *q(m+l) 

g 

= TT^f, v (x n ) (13) 

1=1 

= ?E, q (x n ) 

We note that in the product of the operators in equation (|T3|), the operator with the 
lower value of I comes to the left. The Markov chain defined by the random variable, 
{X n } is ergodic and according to Frustenberg's theorem, a unique invariant measure 
dv E {x) for this process will exist on R. Furthermore, this measure will have density, 
i.e. cIi>e{x) = (f>E(x)dx, when 4>e{x) G Li.(R, dx). In other words belongs to 

the class of nonnegative Lebesgue integrable functions on the compactified real line. 



3 Equations for the invariant measure 4>e(x), the 
Lyapunov exponent and the Integrated density 
of states 

Before defining the equation for determining we define an operator, Vo Te such 

that 

( V °T E f)(x) = 7TT^/(^) (14) 
{hi — xy hi — x 



This is the logical extension of the definition given by Bovier and Klein [12). This 
generalization is useful for the more general problem as it can be seen in the forth 
coming discussion. Some important properties of Vo Te which will be used here are 

(i) v «Te- v v °T E - u = T E - v Te-u 



;ii) v °TE l £ V °T E = ^±x\ including E = 0. 



(iii) y °T y °T = I 
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(iv) VoT E = e~ E £ Fo T , 

(v) ( V °T E 1 f)(x) = ^f(E-^), 

If du E (x) defines the invariant measure with respect to the process described by 
{X n }, introducing equation (13) in equation (10) we obtain 

/ du E {x)f{x) = E / dv E (x)f(P E>q (x)) (15) 

for all bounded measurable functions. As usual E denotes expectation with respect 
to the probability distribution of {e qm } and {p qm }- In our model calculations we take 

P{e) = \s{e) + \6{e - 1) (16) 

Since dv E {x) = cf> E (x)dx, introducing this relation in equation QI5D and making ap- 
propriate change of variables we obtain 

/ <f> E {x)dx = E / dx f[ ?(9 '°^_ e(?j0 (j) E {x) (17) 



where e(q, I) and V(q, I) are defined respectively by equation(|TT]) and (|12|) by setting 
Pqm = P and e qm = e. Again in the product of these operators, the operator with lower 
value of I comes to the left. Since equation ( 17|) should hold good for any arbitrary 



bounded measurable function, /, we obtain after averaging over e and p 

4>e{x) = [\t e + \f[ v ^-^T E _ eq _ t + - n v ^T E ^J<j> E {x) (18) 
z 4 1=1 4 1=1 

As an example we consider the SRTM. For this model q = 3, €2 = eo = v and e\ = u. 
The hopping between nearest-neighbors in the guest cluster is Vq. Since the cluster 
has a inversion symmetry, two products are identical. So, we have 

M*) = + V °T E - V V °T E _ U T E _ v ](j) E (x) 

= ^[T E + T E _ v T (E _ u)/V 2T E _ v ](f) E (x) (19) 

Before computing the complex Lyapunov exponent, 7(E) , we note that for disor- 
dered systems including systems under study, due to the subadditive ergodic theorem 
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D| the limit in equation(6) exists and is independent of the realization of the 



disorder, for all most all realization. In other words, j(E) is self averaging, 

E7(£) = i{E) (20) 

where E denotes ensemble average. For effective calculation of j(E), we first write 

7(E) = lim ^ N (E) (21) 

where N = qM is the total number of sites in the chain. Since 

i 

Zqm+l = Yl ^>E,e(q,l-k),V(q,l-k)( Z Q"i) 
k=l 

= P E ,i{x n ) (22) 
[see equation (13)], from equation (6) we obtain 

^ M I AI 9-1 1 M 9-1 ~ 

*< N W = ^ 111 Xn + ^ ^ 111 ~ E ln II V q n+l,qn+l-l (23) 

i l lvl n =0 n=0 1=1 n=0 1=0 

The last term in equation (^) when averaged over all possible realizations of the 
sample (over e and p) yields 

^ m q-i ^ q-l 

< ~T7 E 111 II Vgn+l^n+l-! > e , p = — In J] VU,Z (24) 
Q M n=0 1=0 Z( l 1=1 

In the limit iV — > oo we apply the subadditive ergodic theorem to equation (p3|) . This 
in turn yields 



7(.E) = - / dv E (x) hxx 
q Jr 



+ -]TE [. dv E (x)\nP E>l (x) - ^InllV^ (25) 

where E denotes expectation over {e} and {p}. Furthermore, introducing dv E (x) = 
4> E (x)dx and making appropriate change of variables in the second integral of equation 
), we obtain 



■y(E) — - [ dx<p E (x)\rix 
q Jr 

+ -£E / cfelnx [J 9{q > l - m) T E -e(q,i- m )Mx) ~ ^^YlV^u (26) 

9 1=1 jR m=0 Z Q 1=1 
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We again point out that in the product of the operators, the operator with the lowest 
value of m comes to the left. Furthermore, the product terminates at m = I. Equation 
(26) when averaged over e and p yields 

If 1 9-1 f 9-1 

j(E) = - I dx4> E (x)\nx + — V /. y^dx\nxT l E (j) E (x) 
q jr 2q y— j JR 

+ J R dx\nx n Yl - mM T E _ n _ m <\> E {x) 



I g-l . i 
H V / dxlnx TT v q -i-i+m, q -i+m Tj 



g-i „ i 

e q _ l+m ^<pE[X) 



(27) 

Before applying equation fl2T|) to the SRTM we first observe that 

/ dx \nx Vo T E _ u T E _ v <fi E (x) — In Vn = / dxlnx T(e-u) T E _ v (j) E (x) (28) 

Decomposing the complex Lyapunov exponent j(E) to real and imaginary parts, we 
obtain for the IDOS, N(E) of the SRTM 

N(E) = 1 - -Ini7(£) 

71 

= 1 - - / dx[2 + T E + Tj + T E _„ + T( B _ U ) T B -«]0B(a;) (29) 

J-oo v 2 



Furthermore, for the Lyapunov exponent 7(-E) we get 
7(E) = Re7(£) 

1 Z" 00 

= - dx\ny(x)[T E _ v T ( _E-u 1 T E _ v -T E ](j) E (x) 

O J-oo vj— 

1 r 00 

+- / dx\n \ x \ [T E - V + T(E-u) T E - v ]</) E (x) (30) 

O J-oo v? 



To obtain the first integral of equation ( |5UD we define y(x) = \J x 2 — Ex + 1 so that 
| a; |= y(x)/y(E — -). This in turn yields 

In I x 1= lny(x) — lny(.E ) (31) 

x 

and the integral decomposes into two integrals. In the second integral we replace x by 
1/(E — x) so that y(E — -) — > y(x). Then we combine these integrals and introduce 
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the governing equation of 4> E (x) [equation (|T9|)1 |TjJ for the final result. We further 
note that the effect of the guest cluster on 'y(E) appears in the second integral of 
equation fl30|) . So, the transformation required for this integral will depend on the 
structure of the guest cluster. For example, consider the example of the RBA. Here, 
v = and u — (1 — V 2 )u. For this case we simply introduce equation ( pT|) in equation 
( PU| ) and after some algebra we obtain 

1 f°° 

j(E) = ~ dx\ny(x)(T(E-ui -T E ) Te4>e(x) (32) 

O J-oo v 1 



equation ( j32|) immediately yields 7(0;) = confirming the result of ref. 0. As a 
second example, we consider u = (1 + V 2 )v in the SRTM. In this case the system has 
an exceptational energy at E = v. To obtain the behavior of ~f(E) in the neighborhood 
of v we transform equation (|30D to 
1 f°° 

7(E) = - dx\ny(x)[T E -vT(E-u)T E _ v -T E ]<f> E (x) 

O J —00 

1 r 00 

-- / dx\ay(x)[l - T E \ T [T^-u + T (E -u) T E _ v ](j) E (x) (33) 

O J— 00 

Since T 2 = I and T T = T^ 1 , from equation ( |3"3"D we obtain 7(f) = which is 
in agreement with the result in ref 0. When Vq = 1, this system also has another 
exceptional energy at E — 2v. The calculation of j(E) for this case is facilitated by 
converting equation (|30|) to 

1 f°° 

l{ E ) = 7 dx m y{x)[T E -v T E -iv T E _ V - T E ](fi E (x) 

O J-00 

1 r 00 

+- / dxlia \ x \ [T E _ V - T Q T E _ 2 „ T E - v ]<p E (x) (34) 

We note that when E = 2v, T E - V T E _ 2 v T E _ V = T v T T v = T 2v and T T E _ 2v T E - V = 
T v . So, j(2v) = 0. This result has also been obtained in ref. H]. There is another 
important motivation behind bringing lny(x) in the integrals. When E = 2cosa7r, 
these integrals can be easily carried out by mapping the compactified real line on the 
circumference of a circle of circumference, tt(Si) through the transformation, 

x = Sin(g + a7r) , (35) 
sin 

when x G R and 9 G Si. As it will be shown later, this transformation will also 
simplify the perturbative calculation of 4> E (x). 
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4 Exceptional energy 



For the purpose of clarity we discuss this aspect in reference to the SRTM. We already 
obtained the governing equation of 4>e(x) for this model [equation ([19])]. We note that 
for the RBA, at E = uj equation ([B^) reduces to 

<f> u {x) = Tl <f> u (x) (36) 

For the SRTM with u = (1 + Vg)v we obtain at E = v 

M^ = \[T^ + T- 1 }M^ (37) 
and for V = 1 at E — 2v, we get 

<p2v(x) = ]^[ T 2v + T v To T v \(j)2v{x) 

= l -[Tl+l]T 2v <j> 2v {x) (38) 
We note now that Te f — f always has the normalized solution 



-E 2 /A 

and f(x) E L\(R,dx) iff | f |< 1. For E = 2 cosa7r, if a is irrational, this solution 
is unique [12|, [L7|| . 

It is important to note that the solution <Pe{x) of equation ([36]), (0) and fl3~8|) 
is f'E s {x) when Eg = uj, v or 2v. Keeping this in mind, we define the exceptional 
energy as the energy at which <Pe s = fE s { x )- This definition has also a physical 
origin. Consider two periodic structures, one from the host cluster and the other 
from the guest cluster. The periodic system from the guest cluster will in principle 
form q bands. Intersection of the DOS of any one of these bands with the DOS of 
the periodic system from the host cluster determines the exceptional energy (Es)- 
Since the invariant measure density of the host system is given by equation (|3~9]), the 
definition ensues. Then, to derive the equation for exceptional energies of the SRTM 
we need to solve 

f E (x) = T E - V T(e-u) T E _ v f E (x) (40) 
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It is relevant at this point to note that if 

9( x ) = ~ 7 h\2, 2 ( 41 ) 

then Tg g{x) is also another Lorentzian distribution centered at bi with a half-width 
ai, where 

b 1 +ia 1 = E- — ^- (42) 
b + ia 

Hence, from equation (|40| ) we obtain 

S = E~v- E _ u 1 i (43) 



where £> = ^ + — E 2 /A. equation fl43|) in turn yields a quadratic equation in S. 



The required equation is obtained from either KeS = ^ or from ImS = J 1 — E 2 /4. 
Both will yield identical equations. This procedure for the SRTM yields 

w(£ 5 - v) 2 - [v(u -v) + (l- V 2 )}(E S - v) + [u - (1 + V >] = (44) 

This equation has already been obtained in ref fTj by considering the reflection co- 
efficient of a single guest cluster embedded in the host lattice. Full analysis of this 
equation can also be found there. We note that for v — and u — (1 — Vq)uj, from 
equation ( [4*4]) we obtain E$ = uj. On the other hand if we take u — (1 + V^ 2 )t> , from 
equation ( |44"| ) we obtain £' Sl = t> and -E S2 = (1 + Vq)v + (1 — Vq)/v. From these 
we can immediately see that two exceptional energies in this case will merge at v if 
Vq = jz^2- For Vo = 1, similar situation can be obtained by setting u = [v+ 2 ^~^J — —} 
and | v |< 1. For this case, two exceptional energies will coincide at Eg = (u + v)/2. 

We further note that the procedure outlined here can be applied to complicated 
one dimensional chains like polyaniline, polythiophene etc. We need to apply first 
the real space renormalization procedure to these systems to obtain effective one 
dimensional chains ||. Then we can apply this procedure for finding exceptional 
energies in these systems. 

When the guest cluster is asymmetric /e(x) has to satisfy two equations for two 
different orientations of the cluster. This in turn will yields two equations for 5* and 
these need to be identical. This develops further constraints in the parameter space. 
Because of this extra constraint on the equation for S, the probability of obtaining 
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exceptional energies with asymmetric guest cluster, particularly if the cluster size is 
relatively large, is negligible. One example of an asymmetric guest cluster having 
an exceptional energy has been worked out in ref. 0. This shows the importance 
of symmetry in the guest cluster for obtaining exceptional energies. This aspect has 
also been discussed in ref. 0. 

5 The perturbative calculation of the invariant mea- 
sure density, </>e(x) around the exceptional en- 
ergy, E s for the SRTM 

The procedure for the calculation is well documented in the literature. For E = 



2 cos air, the calculation is facilitated by transformation given by equation (35) . Since, 
for the invariant measure density we should have 

ct) E {x)dx = h E (9)d9 (45) 

We define an operator J a such that 

(1 T 

(Ja<t> E ){0) = Mx)^ = h E (6) (46) 
We also define another operator, r a = J a T E J" 1 such that 

(r a g)(9) = g(9 - an) (47) 
Furthermore, if g(x) = ^j, we then have 

(J a G)(9) = g(x)^ = ~(J Q f){9) (48) 

where (J a f)(8) = f(x)%. It then follows that 

j d j_ x d d6 d sin 2 9 
dx 01 d9 dx d9 sin air 

We now consider the RBA and put E = to + A, where E = 2cosa7r and a is 
not belong to Q JT7J. It is transparent from equation ( fL9|) that the perturbative 
calculation will be less cumbersome if we replace u> by (E — A) . This yields 

(f> E , e (x) = \[T% + T E T E _ e T E ]4> E ,e(x) (50) 



15 



when e = — A(l — Vq)/Vq. So, the problem effectively reduces to the calculation of 
the invariant measure density of a random binary mixture of two trimeric clusters. 
The central energy of the guest cluster is e and all nearest neighbor hoppings are 
unity. We again emphasize that we are seeking (f>E,e( x ) G L^_(R,dx) for | E |< 2. To 
obtain the perturbative solution of (pE j£ (x) around e = 0, we write 

oo n 

and 4>e\ x ) satisfies $^4$' \x)dx = 8 ni0 . We further note that 4>° E (x) = ^r ^^S ;- 
So, 

,4(9) = = _£2=<HI> (52) 

at) 7r 

when sgn(a7t) is positive or negative depending on whether a is positive or negative. 
We know that according to Anderson's theorem, all eigenstates of the system will be 
exponentially localized for | e |> 0. Furthermore, the Lyapunov exponent, j(E) is 
the inverse localization length of the eigenstate at E. So, physically r y e (E) > for 
E G (—2,2). We also know that ^f e {E) is a continuous function of e particularly if 
| e |<< 1. So, we can expand J es (E) in a Taylor series around e = to obtain 

oo n 

le(E) = £ ry7n(£). (53) 

n=0 n - 

Since, 'jo(E) = 0, j(E) reaches the minimum value at e = 0. So, for | E |< 2, we 
must have 7i(-E l ) = and the leading order term in the Taylor expansion will be 
0(e 2 ) with r )2{E) > 0. We also note in this context that when | E |= 2, although 
7o(-E') = 0, however the analytical continuation of ^(E) for | E (> 2 will yield negative 
^{E). This naturally follows from the constraint on <Pe(x). So, in this case 7i(£") 
will be non-zero. We also see from equation(|3~2"D that ^ m {E) has no contribution from 
4>e { x )- So, to calculate the leading order term of j e (E), the knowledge of 4>e\x) will 
suffice. Similarly for the density of states at Eg, we do not need more than <p^\x) 
[see equation (29)]. From equation (|50|) and equation (^l|) for (J a (f)^)(6) = h^\6) 
we obtain 

H ^3^(1) 1 d sin 2 (g-mr) 3 

The procedure for solving this type of equations is to expand h E (9) in the Fourier 
series in (-^-, f ) and then calculate the coefficients from the governing equation. This 
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procedure when applied to equation ( |54"D yields 

h<&\6) = (zll + (55) 

47r | sina;7r | sin3a7r 

The second case that we consider is Vo = 1 and u — 2v. Here, v and 2v are 
exceptional energies and we have already shown that ^(Eg) = for Eg = v and 2v 
[See equation (33) and equation (34)]. Since Eg = E — A, from equation ([19]) we get 
for | v |< 1, 

4>e(e s ),x(x) = 1[T| + T B _ El T E _s 2 r B _ Sl ]0 E(£ ; s ) jA (x) (56) 

where Si = — A and S 2 = 2(_E — A) if Eg = v. On the other hand, for Eg = 2v, 
we have Si = and S 2 = E — A. So, here the effective site-energies of the guest 
cluster depend on the energy under consideration. When viewed from a broader 
perspective, such a situation does arise in the analysis of complex one dimensional 
chains like polyaniline, polythiophene etc. But, the situation there is a bit more 
complex. In relation to the localization of eigenstates, arguments presented for the 
previous example also hold good here. States which are delocalized for A = will be 
exponentially localized for A 7^ 0. Again, the Lyapunov exponent, "f\(E) is expected 
to be 0(A 2 ). So, to obtain the first order term in the perturbative calculation of 
4>E(E s ),\{x)i we write 

{J a (t>E(E s ),\){0) = h E{Es)>x = -r hP {Es) (9) (57) 

n=0 n - 

Then from equation ([56]) and equation fl5T|) for Eg = v and 2v we get 

(I - 2r Q + r a 4 )h« , (6) = - sin 26 (58) 

v ' 7r I smew I 

and 

q >. . (w / „n i? 2 sin 2$ 

respectively. Since we plan to calculate only the IDOS for these cases, we solve for a 
new function, f E ( Es ){6). For Es = v and 2v 

f E(v) (x) = [l-T E - T E 2 - Te 3 ]^ (x) (60) 
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and 

Now from the governing equation of fiE(E s )(Q) [i-e equation (p8|) and equation(|59l)1 we 



f E (2v)(x) = [2 + T E + T E 2 ]4> ( S 2v) (x) (61) 



find that 

x? E 2 sin2# . . 

I - r a f E(Es ) (0) = ^ , | . 1 62 

v S j\ C(E s )7r I sma7r I 



when C(S) = 1 and —4 for E$ = v and 2v respectively. This equation fl6"2"|)can be 
solved by the standard procedure and we obtain 

fE(E s )( x ) = (Jaf E (E s )){Q) 

:<Pe\ x ) ~ o^/^ ■ -^e 0) 



AC(E S ) sin 2 an * K ' 2C{E S ) | sinmr 

Finally we note that when u = (1 + Vq)v and Vq{1 — v 2 ) = 1, two exceptional 
energies merge at E = v [7(f) = (equation (33)]. For this case the equation for 
h E {y),\{Q) [~to note merging of two Eg 's] is 

(I-2r a + r*)hg^(e) = (64) 

The solution for this equation is hSL(0) = constant. Since 

^ h E{ v)A6)dO = F' 2 htLWe = 1 (65) 

-7r/2 J — tt/2 

by construction, we need 

h%, ) (0) = (66) 
to satisfy the constraint imposed by equation 



6 Calculation of the Lyapunov exponent 

We consider here the RBA and the SRTM with a degenerate exceptional energy. 
We choose the first problem because the guest cluster has different hopping elements 
inside it. On the other hand, the second problem allows us to investigate analytically 
the effect of a degenerate exceptional energy on j(E) and the IDOS. 
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(I) The RBA: 



Since T E _ e = e e dt T E , from equation (32) and equation (55) we obtain 

1 f°° d 
7l (£) = - J_Jx\ny{x)—<t> E {x) = (67) 



and 



1 r°° d 1 r°° d 

12 {E) = -|^^lny(a;)(-) 2 0^) + -y_ oo rfa;ln 2/ (x)-T|0^a;) 

1 rl 2 

= o / d6 fin I sin ait I — In I sin 9 II (cos 2$ — cos 4$) 

6vr sin 2 «7r 7- tt/2 11 ; 



37r | sin 
1 



- r f' 2 d9 [In I sinmr I - In I sin 9 11-^- sin 2 6» r 2 h { E ] (6) 

in cot | J-Tr/2 L 1 1 1 u d9 a E y ' 

(68) 



12(4 - E 2 ) 
So, for this case we have 

when | uj |< 2. 

(II) The SRTM with degenerate exceptional energy : 

We define an operator, 0(E, A) which is the total operator operating on <f> E (v)(x) 
in equation (33). 

0(E, A) = T\T- E +\f(E,\)T\ — T E — (1 - T E )T \T\ - T_ E +\f{E,\)T\\ (70) 
where /(£, X) = 2- E 2 + 2E\ + A 2 . To expand 6(E, A) around A = 0, we write 

oo 

6(£,A) = ^A"0„(£) (71) 

n=0 

Now using the standard procedure, we obtain 

O (E) = (72) 
1 (E) = -\j-{T E + T^ + {2-E 2 )\] (73) 
0,(23) = ^) 2 ^-T- + (2-i, 2 ) 2 I ] -2^ + (2 _ i , 2) ^ 2 



2^/1 - g/4 i 

7r dx dx (p E (x) 



--WT^Te 1 (74) 
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From equation (33) and equation (66) we find that 

4-E 2 r°° , , , , d 



7l (£) = — J dx\ny(x)^° E (x) = (75) 



and 



1 /■oo A 9 roa 

l2 (E) = ~(2-E 2 ) 2 dxlnyWi-nUx) - - dx\ny(x) 
6 J-oo dx 3 J-oo 

(2-e 2 ) r°° d d 2j0 . . 

+ - / dx\ny(x) — — x 2 (f>° E (x) 

J —oo dx dx 



(2-E 2 ) 2 W2 



6-7T sin 2 an J-n/2 



/ d9[\n | sinew | — In | sin 6 1 |] (cos 29 — cos 49) 

J-tt/2 



(2 — E 2 ) r' 2 

H — ^ '— \ d9\\n I sin an I - In I sin 9 II [cos 2(9 + a7r) 

3tt sinew L/-W2 11 I I ul v ) 



(2 - £ 2 ) /-^ 

-7T/2 

-cos(46> + 2cnr)] 

(76) 



Since 72(f) = for this case, in principle further calculation is needed to find the 
leading non-zero term in the expansion of j(E). However, the leading term can be 
obtained through a simple argument. For any arbitrary Vq, the original system will 
yield two exceptional energies, Es 1 and Es 2 provided | Es { |< 2 for i — 1,2. If 
I v I < 1, so that Vq is positive, for E$ 2 to be the second exceptional energy we need 

1 -^Pr<V 2 < 1 -±M (77) 
1+ I v I ~ ~ 1- I v I v ; 

Now from equation (30) we find that 

7 (£&) = -J_Jx\n\x\[T Es2 _ v + T E -l 2 _ v }<P% 2 (x) 

= (78) 

Furthermore, 7(f) = 0. Since 'y(E) reaches the minimum value at v and E$ 2 , it must 
possess a maximum in between these points. As we bring E S2 towards v by tuning 
Vo, this maximum also moves towards v and the value of 'j(E) at the maximum 
simultaneously reduces. In the limit when V^ 2 (l — v 2 ) = 1, two minima and a maxima 
merge at v. So, it is an inflexion point of 7(E) and \e=v should be zero. This is 
precisely obtained. Since for | v |< 2, 7(f) must be a positive semidefinite quantity 
with the leading order term determining the sign, J^r \e=v must be zero and we 
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should have ~f(E) ~ (E — v) 4 . This prediction can be tested by rigorous calculation. 
We shall, however, present an alternative justification. Since around the exceptional 
energies, the system behaves like a weak disordered system, around these energies 
j(E) ~| r(E) | 2 , when | r(E) | 2 is the reflection coefficient of a simple guest cluster 
in the host lattice (ref. [0]). For this system | r(E) | 2 can be found in ref. [J]. The 
Taylor series expansion of | r(E) | 2 around v for | v |< 2 yields 

I r(E) | 2 ~ v2 £~^\ e - v f + 0[(E- v) 5 } (79) 

This is consistent with our arguments. 

Before concluding this section we show that for the SRTM with Vq — 1, u — 2v 
and | v |< 1, 71(f) and 72 (2v) are indeed zero. Consider first the case of v. From 
equation (33) we obtain 

1 f°° , , , , d 



poo H 

li(E) = — \ dx\ny(x) — [Tz 1 + T E + 2]0 

J —00 uX 



2 r°° d 
= ~2 J-oo dxlny ^dx^° E(yX " > 

= (80) 
On the other hand for 2v, from equation (34) we obtain 

7 i(£) = --/ dx\ny(x) — [T^ + T E + 2(--x) 2 T E }cf> E (x) 
12 J-00 ax 2 

+ ^ / dx In I x I ^-x 2 T E / 2 <f E (x) 
o J-00 dx 

E 2 f°° d 1 E 2 r°° d 

= ~24 7-oc ^v(^e(x) - g(l - T ) /_ ^ln I x | T T E/2 <f E {x) 



= (81) 

This is so because both integrals involve odd function of x. Hence, around Eg = v and 
2v, 1(E) ~ -^{E — Eg) 2 . The direct calculation of 72 for these cases albeit possible, 
is however quite complicated. But as mentioned earlier a good estimation of 72 can 
be obtained from the reflection coefficient of the single guest cluster. 

7 Density of states at the exceptional energy 

(I) The RBA: 
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From equation (29) we get 

N(E = u + \) = 1-/° dx(fP E {x)- £ - [ dx(l + T E + T E )^\x) 

J — oo O J — oo 

- £ -<t> E (0) + O(e 2 ) (82) 
Again from the governing equation of h E \6), i.e, equation (54), we obtain 

f E (9) = (I + r a + r 2 a )h E \e) 
cos(2# — an) 



An sin an | sin an \ 

This equation, in turn yields 

df) 

M = fM- 



(83) 



8sin 2 a7r(E 2 - l)^ v ' 4(E 2 - 1) dx 

&\ x ) ( 8 4) 

4 | sinmr | (E 2 — 1) v ; 

Required integrals for the calculation of N(E) can be easily performed. The DOS at 
cu, p{uj) is 



dN(E) 



\E=u 



dE 

= 1 | ^ - ^ m 
nV^^ 6V 2 nV^^ 

(II) The SRTM with a degenerate exceptional energy 

We have already shown in this case <f) E \x) = 0. Furthermore, to obtain the DOS 

at E = v, we need the coefficient of (E — v) — A in the Taylor series expansion of 

N(E). So, we write from equation (29) 

N (E = v + \) = I ~lf° <t> E {x)dx-\ f° dx[T E _ v + T^uT E _ v \(tP E {x) 

12 , E 1 , E-lv 
H tan : H tan 



2 3tt V4-£ 2 6tt V4 - £ 2 



where 

E — u E-2 

b(E) = 



V 2 2{y 2 -Ev + 1) 
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and 

a ^>~ 2(v 2 -Ev + l) 
Now expanding Nq(E) in the Taylor series around E = v , we get for the DOS, p(v) 

( ) 1 ^ (87) 

' 7T\/4 — 67r\/ A — V 2 

We consider now the case of the SRTM where Vq = 1 and k = 2t>, and | t> |< 1. 
For this case the DOS at E = v and 2v has been calculated numerically. We present 
here the analytical results. We write for these cases 

N(E = E s + A) « JV (£ = £ 5 + A) + —Nx(E = E s + A) + 0(A 2 ) (88) 

6 

when iV (i? = Es + A) is obtained from equation (|86|) . 
Furthermore, 



iV^E) = / dxf E{Es) (x) (89) 

J b(E) 

and 6(i?) = £7 and for Es = v and 2v respectively. fE(E s ) * s given by the equation 
fl63|). After performing the required integral and combining the coefficient of A in the 
expansion of N(E) we get 

p(v) = ; H ; H V ; ' (90) 

vrv/4^2 4^743^2 24vr v / 4^2 

and 

1 t» 2 

Finally combining the results of the two sections we find that the number of states 
having localization length superior to the sample size is ~ p{Es) '(Es) ] 1 ^ 2 " wri ere 
M is the size of the sample, n denotes the degeneracy of the exceptional energy. 
Consequently, the mean square displacement of an electron should go as t 2 " 1 with 
7 = (1 — j-). This prediction matches very nicely with exponents obtained from 
numerical simulations ||, [TIJ f2l| . 



8 One dimensional correlated disordered system 
as an effective Lloyd model 

We have already proved that the invariant measure 4>%{x) of the system at the ex- 
ceptional energies is a Lorentzian distribution centered at Es/2 with a half width 
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^/l — Eg/4. We have also seen through examples that the DOS at Eg in the near 
perfect limit to a good approximation is the DOS of the perfect system. It is also a 
well established result that the Lyapunov exponent, 7(E) around Eg can be approx- 
imated to a fair degree by | r(E) \ 2 where | r(E) \ 2 is the reflection coefficient of a 
single guest cluster in the host lattice and | r(Eg) \ 2 = 0. 



The Lloyd model Eq] on the other hand is the uncorrelated site disordered An- 
derson model where the probability distribution of the site energies {e n }, P(e n ) is 

p (en) = - ? (92) 

The invariant measure 4>e{x) for this model is 

1 e 
n(x-e* ) 2 + et 



b(x) = ~ - A ■ -2 (93) 



where 



4 = f + ^{E + ie^-A = | [1 + A] (94) 



e! = ^— + [ei + VI] (95) 



and 



a = i±jLbto^±ig^ (96) 

6 = 73+P (97) 

We note that all nearest neighbor hopping matrices have been assumed to be unity. 
The Lyapunov exponent and the DOS for this model can be found in the literature 
We simply quote the results. 

4cosh 7 (e) =| 2 + E | [1 + ^f^+ \2-E | [1 + j^f ?x] 1/2 (98) 

and 

Now in the limit e x -> 0, we find that e*(£) = f , ^jg! = ±, eJ(S) = - E 2 and 
~^nP~ = ~ 2Vi-E i ' Consequently, we obtain 



1 Jl-£ 2 /4 
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l{E) = 



(101) 



and 



p(E) 



1 



(102) 



ny/A - E 2 



So, all the characteristic features of one dimensional correlated disordered systems 
around the exceptional energies are recovered. To develop an effective Lloyd model 
for one dimensional correlated disordered systems, we need then l\ = f(E) such that 
f(Es) = for | Eg |< 2 and f(E) should also contain in it the information about 
the guest cluster. We further note that for Z\ ~ 0, equation ( pgp yields ■y(E) ~| 1\ \. 
But we have shown that for a nondegenerate exceptional energy, ^(E) ~ (E — E$) 2 . 
So, f{E) should also satisfy this condition. In as much as | r(E) | 2 satisfies all these 
criteria, we propose that l\ =| r(E) | 2 . This proposal in turn implies that 



1 V4 - E 2 


r 


(E)\ 


2 


7T e 2 + (4 - E 2 ) 


r(E) 


1 



(103) 



in the effective Lloyd model. Finally we note that many methods are developed 
to study the electrical conductivity of uncorrelated site disordered system. So, these 
methods can be applied to the systems considered here through the proposed map- 
ping. 

9 Summary 

The behavior of electronic states of one dimensional correlated disordered systems 
around exceptional energies (E$) is studied analytically using the invariant measure 
method. The RDM is the simplest example in this category and this has been studied 
by this method by Bovier. The basic approach of Bovier is generalized thoroughly 
and rigorously to take into consideration more structure in the guest cluster. The for- 
malism is further elaborated by applying to the SRTM. Another useful contribution is 
the alternative mathematical definition of the exceptional energy from the invariant 
measure. This definition is further substantiated by physical arguments. Further- 
more, from our definition of exceptional energy we obtain an equation constraining 
the parameters of the guest cluster. The same equation has also been obtained by set- 
ting | r(E) | 2 = 0. This clearly shows an intimate relationship between the invariant 
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measure and | r(E) | 2 . This relationship is further highlighted here by mapping these 
systems to an effective Lloyd model. This equation further shows how the structure 
of the guest cluster can be modulated to tune positions of exceptional energies. We 
also obtain through it the condition for a degenerate exceptional energy. Hence, the 
importance of the method is further illustrated. 

In relation to the localization of eigenstate we find expectedly ■y(E) ~ (E—Es) 2 for 
cases with a non degenerate exceptional energy. In the case of a degenerate resonance, 
72 (-Es) along with "fi(E s ) are shown to be identically zero. These results are further 
substantiated by rigorous analytical arguments. Further analytical arguments are 
presented to show that i(E) ~ (E — Es) 4 for this case. Although a system containing 
a degenerate exceptional energy is studied previously by us, this is, however, the most 
rigorous analysis. 

In the RDM p(v) is the DOS of the perfect system at E — v. Further structure 
in the guest cluster is found to be manifested in p(E$) through correction terms. 
However, the universality in the mathematical expression of Ni(E) should not be 
overlooked. This exemplifies further the universality of one dimensional correlated 
disordered systems around exceptional energies. We finally add that the real signifi- 
cance of our work along with that of Bovier on the RDM is that this firmly establish 
the anomalous behavior of one dimensional correlated disordered systems around ex- 
ceptional energies. 
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